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Strongly interacting bosons in a two-dimensional rotating square lattice are investigated via a 
modified Bose-Hubbard Hamiltonian. Such a system corresponds to a rotating lattice potential 
imprinted on a trapped Bose-Einstein condensate. Second-order quantum phase transitions between 
states of different symmetries are observed at discrete rotation rates. For the square lattice we study, 
there are four possible ground-state symmetries. 

PACS numbers: 



One of the most exciting recent advances in the field 
of ultracold quantum gases has been the realization of 
quantum phase transitions in Bose-Einstein condensates 
(BEC's) trapped in a lattice 0, in particular the super- 
fluid to Mott insulator transition |2j. Quantum phase 
transitions typically occur when a topological change 
happens in the ground state symmetry as a function of 
a Hamiltonian parameter at zero temperature They 
play an important role in many models, including the 
quantum Ising model, quantum rotors, and the Bose- 
Hubbard model . The latter very accurately describes 
experimental observations of the superfluid-Mott transi- 
tion in BEC's 

A second fruitful development in BEC's has been the 
study of rotating systems in harmonic traps, in which an 
Abrikosov lattice of quantized vortices appears for suf- 
ficient rotation 0. Rotating BEC's have also been pre- 
dicted to realize an analog of the fractional quantum Hall 
effect 0- More generally, the density, temperature, lat- 
tice structure, strength and symmetry of interactions in 
ultracold quantum gases are precisely and dynamically 
controllable in experiments, with no impurities or dis- 
order. Thus these systems are accurately described by 
relatively simple model Hamiltonians. In this sense they 
are a playground for quantum many body theory. 

In this Letter, we pull together these two cutting edge 
areas of study in BEC's, i.e., rotating systems and quan- 
tum phase transitions in lattices. Specifically, we derive 
and solve a modified Bose-Hubbard Hamiltonian for ro- 
tating BEC's trapped in a two-dimensional (2D) square 
lattice potential with variable filling factor. In addition 
to the superfluid-Mott transition, which can be driven by 
hopping and/or rotation, we find that a quantum phase 
transition occurs each time the rotational symmetry of 
the ground state changes. The possible discrete values of 
the rotational symmetry are determined by the geometry 
of the underlying lattice. Fig.^is a sketch of one way to 
make a 2D rotating lattice.as has already been experi- 
mentally realized at JILA Some theoretical aspects 
of this problem for large filling factors, such as vortex 
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FIG. 1: (color online) Sketch of a rotating lattice. A laser 
beam passes through a rotating mask. A lens focuses the 
resulting lattice beam profile onto a stationary, trapped BEC. 
The lattice potential is imprinted via off-resonant interactions 
between the laser and the atoms. 

pin ning |||, structural phase transitions of vortex mat- 
ter |l(J| . and the single vortex problem have been 
studied previously. In contrast, we will consider filling 
factors of less than unity. In particular, we treat small, 
finite systems where exact solution via diagonalization 
of the Hamiltonian in a truncated Hilbcrt space is in- 
sightful. Recall that, in most experiments, condensates 
typically are comprised of 10 6 atoms, 3D optical lat- 
tices consist of 100 3 sites spaced at 300-500 nm pj, and 
exact solutions for systems of this size are not tractable 
numerically. Our approach builds up a scalable picture of 
the lattice physics from the microscopic interaction and 
lattice structure. 

Consider bosons interacting via an effective 2D two- 
body contact potential of strength g in a 2D lattice po- 
tential V lat (x,y) = V lat (x + jd,y + kd) which rotates 
about the z axis, where j, k are integers and d is the 
lattice constant. In the rotating frame, the Hamiltonian 
is 



H = Jd 2 r^ 

L z = - ih(xd y - yd x ) 
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with rotation frequency f2 and atomic mass M. The field 
operator, y), obeys the usual bosonic commutation 
relations. Using a Wannier basis, Wi(x,y), the field op- 
erator can be expanded in terms of bosonic operators a,-, 

Hx,y)=J2 i a i W i (x,y). (3) 

The usual single-band Bose-Hubbard model 0, 13 is ob- 
tained via the tight binding and lowest band approxima- 
tions The rotational part H L = J d 2 r&Q,L z <& of 
Eq. 0} becomes 

^ = -M/<i J rE<y> 

a]ai(Wj(x,y)(xd y - yd x )Wi{x,y)) + h.c. 

= ihtt E<i,j> Ki^aia] - a\a 3 ) , (4) 

where <i,j > indicates a sum over nearest neighbors. 
The general form of ify is 

ify = I3{nrj/d 2 ) sin ay , (5) 

where denotes the distance from the axis of rotation 
to the i" 1 site, a%j is the angle subtended by the i th and 
the j th sites with respect to the axis of rotation, and 
(3 is a dimensionless constant characterizing the lattice 
geometry and depth. 

The rotating Bosc-Hubbard Hamiltonian in the canon- 
ical ensemble is 

H = -tY i<i j > {a i a) + a\a j ) 

+ |^Ei"i(ni-l), (6) 

where hi = aja^. The hopping energy t and on-site inter- 
action energy U have been explicitly calculated from g. 
M, etc. elsewhere Equation (JSJ) is the usual Bose- 
Hubbard Hamiltonian with the addition of the second 
term, which describes the rotation. 

The expectation value of the current is an important 
ground state observable for rotating quantum systems. 
In the Hciscnbcrg picture the current is given by 

(4) = (i/hd)^^]) (7) 
= hd" " 3 ~ i<lj ' d~ ^ j + lCl] ' 1 

where H+j is the Hamiltonian for sites i,j alone. In 
Eq. |J7J|, the first term is due to hopping while the sec- 
ond is due to rotation. The current is conserved at each 
site i in the rotating frame since the sum of (Jy) over all 
nearest neighbors j is zero. Another useful observable is 
the total current on the lattice boundary C, 

A = (fa/£V)£ ( ij) et7 <4->. (8) 
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FIG. 2: (color online) The total dimensionless current A 
around a lattice unit cell is shown as a function of the av- 
erage total number of atoms n in the system for M1K > t. 
The curves from top to bottom are for increasing on-site in- 
teraction strength U in units of the lattice recoil energy E r . 
For weak interactions, the atoms behave independently; for 
strong interactions, there is a particle-hole symmetry, as evi- 
dent in the bottommost curve, where n = 1 and n = 3 have 
the same A. The connecting lines here are a guide to the eye. 

where we have scaled away the units via a "recoil" en- 
ergy E r = h 2 /Md 2 and all sums over C arc taken with 
the same sign convention as the helicity of f2. We also 
define two number-related observables: n = Et(^*)) 
the average total number of atoms in the system, and 
v = EifX 7 ^ 2 ) — i^i) 2 )/ Ei("i)j the normalized variance. 
Recall that v = 1 for a coherent state, v > 1 for a phase- 
squeezed state, v < 1 for a number-squeezed state, and 
v = for a single Fock state. 

The ground state wave function must be consistent 
with the discrete rotational symmetry of the lattice. This 
follows in general from the point symmetry group Cj 
of the lattice, which describes the rotation of a j-fold 
symmetric object, e.g., a j-sided regular polygon. The 
number of possible rotational symmetries is equal to j, 
e.g., for the square lattice we consider j = 4. The quan- 
tum mechanical rotation operator about an angle 27T / j, 
7Z(2Tv/j), then commutes with the lattice potential part 
of the Hamiltonian. This means that all non-degenerate 
ground-states are simultaneous eigenstates of energy and 
of the discrete rotational symmetry operator. Further- 
more, [R,(2tt / j)Y , rotates the system through an angle 
2tt and therefore must be the identity operator for the 
wave function to be single- valued. The eigenvalues of 
1Z(2tt/ j) can thus only take the values e 12 *™ 1 ^ where m 
is an integer and m € {0,7 — 1}. In the case where 
the total particle number is commensurate with j, the 
eigenvalues of [7?. (2 71-/7')] can be degenerate, leading to 
interesting effects which will be discussed later. 

For simplicity, we will focus our discussion on the 
strongly interacting case, achieved experimentally either 
by a Feshbach resonance or by turning up the lattice po- 
tential height. In this case one can make a two-state ap- 
proximation which prevents multiple occupancy of 
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FIG. 3: (color online) (a) index m of the eigenvalue e l27rm / 4 
of the discrete rotational symmetry operator 1Z as a function 
of rotation rate, showing the discontinuous changes in the 
ground state symmetry. The rotation fl is in units of E r /h, 
with total number of atoms n = 1, hopping t/E r = 1, and 
/3 = 4.93 1 (as appropriate for a sinusoidal square lattice), (b) 
The derivative of the ground state energy with respect to Q 
in units of h, showing discrete jumps at the points of level 
crossings. 

the same site. Such a Fock space eliminates the inter- 
action term U in Eq. However, there are effective 
strong interactions due to atoms being unable to cross 
each other; i.e., one has a system of hard core bosons. 

Having defined the observables, we first consider the 
case of a single square unit cell, i.e., a 2 x 2-site lattice. 
In order to assess the two-state approximation, we allow 
Fock states with up to three atoms per site: however, 
we keep the average number of particles per site at unity 
or below. We exactly diagonalize the Hamiltonian and 
find the ground state. Due to the competition between 
the hopping and rotational energy terms in Eq. (jBJ), ro- 
tation affects the system only when MIK > t. For a 
total number of atoms n <E {1, 2, 3}, m = for MIK < t 
and m = 1 for MIK > t. As shown in Fig. for non- 
interacting atoms (U = 0) this corresponds to A = 
and A = 2n, respectively. As the interatomic repulsion 
is increased, A decreases to non-integer values. Although 
m is quantized, A is not. However, for t/U <C 1 and 
Hfl/U <C 1, i.e., for very strong interactions, or a very 
strong lattice and small rotation, the allowed values of 
A return to those given by the two-state approximation. 
The two-state approximation (zero or one atom per site) 
is adequate to study strongly interacting systems. 

Consider next a 4 x 4-site square 2D lattice, which 
consists of 9 unit cells, in the two-state approximation. 
The numerical study shows three main results. (1) Sec- 
ond order quantum phase transitions occur each time 
the symmetry m of the ground state changes. (2) As 
the rotation rate increases, eventually the maximum cur- 
rent is achieved corresponding to a phase difference of 




FIG. 4: (color online) The red arrows indicate direction of 
current flow in the rotating frame for t/E r = 1 and n = 1 
at increasing rotations on a 4 x 4 lattice, corresponding to 
m = 0, 1, 2, and 3. (a) Kl/E r = 0.1; (b) h{l/E r = 0.2; (c) 

nn/Er = 0.4; (d) nn/Er = 0.8. 

7r/2 between the wavefunction on adjacent sites along 
the perimeter. (3) For higher fillings, both the hopping 
t and the rotation fl can drive the system through the 
extensively studied Mott-insulator/supcrfluid transition 
QSCJ01 (not discussed further here). 

First, consider the simplified case of one atom in the 
system. To illustrate the QPTs, in Fig. [3] is shown the 
eigenvalue index m and the derivative of the total energy 
E = (H) with respect to il. Exact energy level crossings 
are observed as indicated by the discrete jumps in the 
energy derivative, corresponding to each transition. The 
corresponding circulation patterns are shown in Fig. 21 
In Fig.^Ja), m = and rotation has not yet entered the 
system. In the rotating frame, the current seems to be 
flowing backwards, i.e., clockwise. In Fig. QJb), m = 1 
and a single vortex enters the system and rests at the 
center. Fig.^c) shows a similar pattern although m = 2 
in this case. Fig. 0fd) is analogous to a tightly packed 
vortex lattice with m = 3; the current directions in the 
center with opposite sign of circulation to the flow around 
the perimeter is a trivial result of packing four vortices 
together on a 2D square lattice with an antivortex at the 
center. 

Finally, these results can be extended for general fill- 
ing. The main difference from the single-atom case is the 
presence of effective strong interactions due to the use of 
the two-state approximation. In Fig. we illustrate the 
number density n^, normalized variance u, eigenvalue of 
discrete rotational symmetry m, and the normalized cur- 
rent on the boundary A/nc, all for a total of five parti- 
cles in the 4x4 square lattice. Note that up, the number 
density along the perimeter, depends on fl. In Fig.J^a)- 
(b) the site-dependence of rid and v are illustrated for 
h£l = t = E r . In Fig. |S{c) one observes that A depends 
linearly on fi; this is qualitatively a similar result to that 
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FIG. 5: (color online) Observables for a 4 x 4 lattice with five 
particles and t/E r = 1. (a) Site-dependent number density n 
for hQ/Er — 0.1. There is noticeable depletion on the inner 
sites, (b) Site-dependent normalized variance v for hQ,/E r = 
0.1. Note that v < 1, so this is a number-squeezed ground 
state, (c) The normalized, scaled current Ac/nc (solid curve) 
and m (dashed curve), all on the boundary as a function of 
S7(in units of E r /h). Between each quantum phase transition, 
where the symmetry of the ground state changes abruptly 
as indicated by the jump in the value of m, the normalized 
current depends linearly on the rotation Q. (d) is the same 
as (c) but describes four particles. 

of one atom. However, max.(Ac/nc) ^ 2, due to interac- 
tions. The other essential features of general filling, i.e., 
the energy level crossings as the ground state symmetry 
changes abruptly, are qualitatively the same as that for 
one atom. For many particles this is a non-trivial result 
indicative of quantum phase transitions. 

Fig. El contrasts results for five particles with those for 
four particles. This situation is qualitatively different as 
there are no jumps in the energy derivative, or changes 
in the symmetry of the ground state, fnstead of exact 
level crossings, avoided crossings are observed. This is 
reflected in the smoothed transitions in Ac/nc- The 
reason is straightforward and illuminating, since for four 
particles in a 4 x 4 lattice the application of TZ(2n/4) one 
time leaves the wave function unchanged. In this case, 
the ground state symmetry must be four-fold degenerate. 
Since the ground state and excited states have the same 



symmetry they can mix and curve crossings do not oc- 
cur. For 2-particles in a 4 x 4 lattice, the eigenvalues of 
72.(27r/4) are instead two- fold degenerate, but since the 
symmetry always switches at curve crossing points, this 
exhibits quantum phase transitions qualitatively like the 
five particle case discussed previously. 

Based on symmetry considerations, we identify tran- 
sitions between rotational ground states as second order 
quantum phase transitions. As f2 increases, the state of 
the system changes continuously while the symmetry of 
the ground state changes discontinuously. This defines 
a second order phase transition 0|. Although we have 
considered small quantum systems, the number of sym- 
metry eigenvalues j = 4 remains true for a square lattice 
of arbitrary size. Identical rotational constraints point 
to the same symmetry values for the rotating 3D lattice 
although the exact nature of transitions will depend on 
dimensionality. 

We note that the appearance of a vortex lattice in the 
discrete system should be observable by expansion and 
interference with a non-rotating system, as in the con- 
tinuous case 0- A second observable has been provided 
by our calculations: the current on the boundary jumps 
discontinuously with increasing f2. The maximal value 
of A at eigenvalue m = 3, indicates a current pattern 
reminiscent of a tightly packed vortex lattice. 
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